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SUMMARY 


This report presents the development of lateral-directional flying qualities guidelines 
with application to eigenspace (eigenstructure) assignment methods. These guidelines will 
assist designers in choosing eigenvectors to achieve desired closed-loop flying qualities or 
performing trade-offs between flying qualities and other important design requirements, 
such as achieving realizable gain magnitudes or desired system robustness. This has been 
accomplished by developing relationships between the system’s eigenvectors and the roll 
rate and sideslip transfer functions. Using these relationships, along with constraints 
imposed by system dynamics, key eigenvector elements are identified and guidelines for 
choosing values of these elements to yield desirable flying qualities have been developed. 
Two guidelines are developed - one for low roll-to-sideslip ratio and one for moderate-to- 
high roll-to-sideslip ratio. These flying qualities guidelines are based upon the Military 
Standard lateral-directional coupling criteria for high performance aircraft - the roll rate 
oscillation criteria and the sideslip excursion criteria. Example guidelines are generated for 
a moderate-to-large, an intermediate, and low value of roll-to-sideslip ratio. 


1.0 INTRODUCTION 

The Direct Eigenspace Assignment (DEA) method (Davidson and Schmidt 1986) is 
currently being used to design lateral-directional control laws for NASA's High Angle-of- 
Attack Research Vehicle (HARV) (Davidson et al. 1992). This method allows designers to 
shape the closed-loop response by judicious choice of desired eigenvalues and eigenvectors. 
During this design effort DEA has been demonstrated to be a useful technique for aircraft 
control design. The control laws developed using this method have demonstrated good 
performance, robustness, and flying qualities during both piloted simulation and flight 
testing (Murphy et al. 1994). 

During the control law design effort, two limitations of this method became apparent. 
First, when using DEA the designer has no direct control over augmentation gain 
magnitudes. Often it is not clear how to adjust the desired eigenspace in order to reduce 
individual undesirable gain magnitudes. Second, although considerable guidance is 
available for choosing desired eigenvalues (Military Standard, time constants, frequency, 
and damping specifications), little guidance is available for choosing desired system 
eigenvectors. Design guidance is needed on how to select closed-loop lateral-directional 
eigenvectors to achieve desired flying qualities. 

The first limitation was addressed by the development of Gain Weighted Eigenspace 
Assignment (GWEA) (Davidson and Andrisani 1994). The GWEA method allows a 
designer to place eigenvalues at desired locations and trade-off the achievement of desired 
eigenvectors versus feedback gain magnitudes. This report addresses the second limitation 
by presenting the development of lateral-directional flying qualities guidelines with 
application to eigenspace assignment methods. These guidelines will assist designers in 
choosing eigenvectors to achieve desired closed-loop flying qualities or performing trade- 
offs between flying qualities and other important design requirements, such as achieving 
realizable gain magnitudes or desired system robustness. 

This report is organized into four sections. A review of lateral/directional dynamics, 
background information on how eigenvalues and eigenvectors influence a system's dynamic 
response, a review of the Direct Eigenspace Assignment methodology, and an overview of 
existing lateral/directional flying qualities criteria is presented in the following section. The 
development of the lateral-directional eigenvector flying qualities guidelines are presented in 
the third section. Concluding remarks are given in the final section. 



2.0 BACKGROUND 


This section presents a review of lateral/directional dynamics, background information 
on how eigenvalues and eigenvectors influence a system’s dynamic response, a review of the 
Direct Eigenspace Assignment methodology, and an overview of existing lateral/directional 
flying qualities criteria. 


Lateral-Directional Dynamics 

The linearized rigid body lateral-directional equations of motion for a steady, straight, 
and level flight condition, referenced to stability axes, are (McRuer et al. 1973) 
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or in state space form 
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where 


P = sideslip angle 

p = stability axis roll rate 

r - stability axis yaw rate 

</> = bank angle 

S a u = aileron control input 

8 ru d = rudder control input 


(2.1b) 


and the prime denotes the inclusion of the inertia terms. As can be seen, the lateral (p) and 
directional (ft and r ) responses are coupled. The primary lateral-directional coupling 
derivatives are: roll moment due to sideslip angle Lp , roll moment due to yaw rate L r , yaw 
moment due to roll rate N p , and yaw moment due to lateral controls Ng . A brief review of 
the physical basis of these derivatives is given in the Appendix. 
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The characteristic equation for this system is 
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There are three classical lateral-directional eigenvalues: a lightly damped oscillatory pole 
referred to as the Dutch roll pole (A^ r ), a first order pole with a long time constant referred 
to as the spiral pole ( X spr i ), and a first order pole with a relatively short time constant 
referred to as the roll pole (A ro //). The characteristic equation can be written in terms of 
these eigenvalues as 

A(s) = k^(s— ^- S p r l )(s — + 2 Cdr®dr s ®dr ) 

v ' (2.3) 
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where k^=Y^-\, ^dr = ~ (0 drCdr + j Cl} dry^~^dr 2 and ^dr denotes the complex 

conjugate of A dr . Approximations for the system eigenvalues in terms of stability and 
control derivatives (McRuer et al. 1973) are given by: 
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The primary lateral-directional control task is control of bank angle with lateral stick. 
The following relationships are developed for lateral stick controlling aileron deflection (<5j,* 
= Sail) with zero rudder input. In the following, the sub-subscript “ail” on the control 
derivatives has been dropped to simplify the notation. The bank angle-to-lateral stick 
transfer function is given by 
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This transfer function can be written in pole-zero form as 
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The following relationships can be written from (2.8) and (2.9) 
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By making the following assumptions (reasonable for most configurations (McRuer et al 
1973)) 


Y r = 0,Yp=0,Y $ =0 


(2.12a) 
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equations (2.10) and (2.1 1) reduce to 


co$=N p 1- 


(2.13) 
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Since p=s(p , the roll rate-to-lateral stick transfer function can be written 
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The steady-state roll rate for a unit step lateral stick input (assuming the spiral pole is 
approximately at the origin) is given by 


• f -i Y a <t> 

Pss Lg 

\*roll )\®dr 


The sideslip-to-lateral stick transfer function is given by 
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Making the assumptions of (2.12a), and that the spiral pole is close to the origin, and that 


Yp +a o =0 


(2.18a) 
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The yaw rate-to-lateral stick transfer function is given by 
S s(k = A(j){ ~ Ns ^ Y p ~ + N§(L p (Yp - 1) - Yp ) - LgN p {Yp - 1) js 2 

+[-Y§(NpL p - LpN p ) + N s (YpL p - Lp(Y p +a 0 )) 

- Lg(YpN p -Np(Y p +ao))} s 
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( 2 . 20 ) 


Making the asumptions of (2. 12a) and (2. 18a), this transfer function can be written as 
r -1 
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( 2 . 21 ) 


Eigenvalues, Eigenvectors, and System Dynamic Response 

The eigenvalues and eigenvectors of a system are related to its dynamic response in the 
following way. Given the observable and controllable linear time-invariant system 


and output equation 


x = Ax + Bu 


y = Cx 

where x e R", u e R m , and ye R l . 

The Laplace transform of equation (2.22a) is given by 

$x(j) - x(0) = A*(.y) + Bu{s ) 


(2.22a) 

(2.22b) 


(2.23a) 


x(s) = [j/„ - A] 1 x(0) + [sl n - A] l Bu(s) (2.23b) 

Solution of equation (2.22a) is given by taking the inverse Laplace Transform of equation 
(2.23b) 

40 = ^“ 1 {[5/„-Ar 1 }jr(0) + ^~ 1 |[j/„-Ar 1 fiu(s)} (2.24) 
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the solution of (2.24) is (Brogan 1974) 

t 

x (t) = e A, x(0 ) + J e M '- z) Bu(T)dT 
o 

and system outputs are 

t 

y(t) = Ce A 'x( 0) + j Ce A(, ~ r) Bu(T)dt 
0 


The system dynamic matrix, A , can be represented by 


(2.26) 


(2.27) 


A = VAV~' = VAL (2.28) 

where V is a matrix of system eigenvectors, L is the inverse eigenvector matrix, and A is a 
diagonal matrix of system eigenvalues. Given this result, e^'can be expressed by 

e* = Ve A, L = £ Vj e Xj ' /,. (2.29) 

7=1 

where Xj is the j^ 1 system eigenvalue, v- is the j^ 1 column of V ( eigenvector of A ), and lj 
is the j^ 1 row of L ( j^ 1 left eigenvector of A ). Equation (2.27) can then be expressed as 

y(r)=X Cv / ><0 ^(°) + Z Cv J fe^ T) ljBu(T)dT (2.30) 

7=1 7 = 1 0 

Noting that 

m 

Bu(t) = ^b k u k (t) (2.31) 

* = 1 

where b k is the k^ 1 column of B and u k is the system input, the system outputs due to 
initial conditions and input u k is given by 


y(r)=XCv y ^ (o /^(0) + XX Cv M je lj( '~ T) u k (T)dr (2.32) 

7 = 1 7=1 *=l 0 

The i^ 1 system output is given by 

m V J x(0) + X £ c,v/A I) dl (2.33) 


7=1 


7=1 *=1 


where c ; is the i 1 * 1 row of C. In the case of initial conditions equal to zero, the i 1 * 1 output is 
given by 


n m t , , , 
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where R .= c v l b. . In this expression, R- , . is the modal residue for output i, associated 
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(2.35) 


Given an impulsive input in the k^ 1 input, equation (2.34) reduces to 


j = 1 k = 1 

and for a step input in the k^ 1 input, equation (2.34) reduces to (for Xj * 0) 
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As these expressions show, a system's dynamics are dependent on both its eigenvalues and 
its eigenvectors. The eigenvalues determine the time constant or frequency and damping of 
each mode. The eigenvectors determine the residues. The residues determine how much 
each mode of the system contributes to a given output. 


For example, for the lateral-directional system given by equation (2.1), time responses 
for a unit step lateral stick input (and zero pedal input) can be written in terms of system 
eigenvalues and residues as (because there is only one input, the third subscript on the R ’s 
has been omitted) 
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and 


vrr / R p<dr 

p A, ’ 




^ = Z-^; 'Fr = Z- Rr ' dr 


x dr ‘ /l 4r A^ r 

where Ixl denotes the magnitude of x and Zx denotes the phase angle of x. 


(2.39) 
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Direct Eigenspace Assignment Methodology 

One possible approach to the aircraft control synthesis problem would be to synthesize 
a control system that would control both the eigenvalue locations and the residue 
magnitudes. Since the residues are a function of the system's eigenvectors this naturally 
leads to a control synthesis technique that involves achieving some desired eigenspace in the 
closed-loop system (eigenspace (eigenstructure) assignment) (Moore 1976; Srinathkumar 
1978; Cunningham 1980; Andry 1983). An eigenspace assignment method currently being 
used to design control laws for NASA’s High Angle-of-attack Research Vehicle (HARV) is 
Direct Eigenspace Assignment (DEA) (Davidson et al. 1992; Murphy et al. 1994). DEA is 
a control synthesis technique for directly determining measurement feedback control gains 
that will yield an achievable eigenspace in the closed-loop system. For a system that is 
observable and controllable and has n states, m controls, and / measurements; DEA will 
determine a gain matrix that will place / eigenvalues to desired locations and m elements of 
their associated eigenvectors to desired values +. If it is desired to place more than m 
elements of the associated l eigenvectors, DEA yields eigenvectors in the closed-loop 
system that are as close as possible in a least squares sense to desired eigenvectors. A more 
detailed development can be found in Davidson and Schmidt, 1986. 


Direct Eigenspace Assignment Formulation 
Given the observable, controllable system 

x = Ax + Bu (2.40a) 

where x e R n and u e R m , with system measurements given by 

z = Mx + Nu (2.40b) 

where z e R*. 

The total control input is the sum of the augmentation input u c and pilot’s input u p 


u = u p + u c (2.41) 

The measurement feedback control law is 

u c = Gz (2.42) 

Solving for u as a function of the system states and pilot's input yields 

« = [/„- GNT' GMx + [I m - GNT' u p (2.43) 

The system augmented with the control law is given by 

x = (A + B[I m - GAT 1 GM)x + B[I m - GN)~ l u p (2.44) 

The spectral decomposition of the closed-loop system is given by 

(A + B[I m - GN]~ l GM)v, = A,v, (2.45) 


t This assumes / > m. For a general statement and proof of this property the reader is referred to 
Srinathkumar 1978. 
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for i = 1 ,...,« where A,- is the i^ 1 system eigenvalue and v, is the associated i*h system 
eigenvector. Let w, be defined by 


w, = [l m — GN] 'GMv i (2.46) 

Substituting this result into equation (2.45) and solving for v, yields 

v i =[X i I n -AT'B Wl (2.47) 

This equation describes the achievable i™ eigenvector of the closed-loop system as a 
function of the eigenvalue A, and w, . By examining this equation, one can see that the 
number of control variables (m) determines the dimension of the subspace in which the 
achievable eigenvectors must reside. 


Values of w, that yield an achievable eigenspace that is as close as possible in a least 
squares sense to a desired eigenspace can be determined by defining a cost function 
associated with the i^ 1 mode of the system 

(2.48) 

for i = 1 ,...,/ where v aj is the i“* achievable eigenvector associated with eigenvalue A, , v e j i is 

the i^ 1 desired eigenvector, and Qd i is an n-by-n symmetric positive semi-definite weighting 

matrix on eigenvector elements f . This cost function represents the error between the 
achievable eigenvector and the desired eigenvector weighted by the matrix Qd i . 

Values of w t that minimize are determined by substituting (2.47) into the cost 

function for v fl( . , taking the gradient of 7, with respect to w t , setting this result equal to zero, 
and solving for w; . This yields 

"i = i A di H QdM^ A di H Qd, v d t (2.49) 

where 

A di = ttd i In~ A r l B (2.50) 

and X tii is the i^ 1 desired eigenvalue of the closed-loop system. Note in this development 
Arf 4 . cannot belong to the spectmm of A. 

By concatenating the individual w,’s column-wise to form W and v a/ 's column-wise to 
form V a , equation (2.46) can be expressed in matrix form by 

W = {I m -GNy'GMV a (2.51) 


^ Superscript H denotes complex conjugate transpose (Strang 1980). 
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From (2.51), the feedback gain matrix that yields the desired closed-loop eigenvalues 
and achievable eigenvectors is given by (for independent achievable eigenvectors) 

G = W[MV a + AW]" 1 (2.52) 


Design Algorithm 

A feedback gain matrix that yields a desired closed-loop eigenspace is determined in the 
following way: 

1) Select desired eigenvalues X ( j i , desired eigenvectors , and desired eigenvector 
weighting matrices Qd l . 

2) Calculate w, 's using equation (2.49) and concatenate these column-wise to form W. 

3) Calculate achievable eigenvectors v a . 's using equation (2.47) and concatenate these 
column-wise to form V a . 

4) The feedback gain matrix G is then calculated using equation (2.52). 


Existing Lateral-Directional Flying Qualities Criteria and Eigenspace Assignment 

A key goal of piloted aircraft control law design is to achieve desirable flying qualities in 
the closed-loop system. A primary source for flying qualities design criteria for high 
performance aircraft is the Military Standard 1797 A - Flying Qualities of Piloted Aircraft 
(and earlier versions - Military Standard 1797 and Military Specification 8785). Using 
eigenspace assignment methods the designer specifies the desired closed-loop dynamics in 
the form of desired eigenvalues and eigenvectors. The Military Standard provides 
considerable guidance for choosing lateral-directional eigenvalues to yield desired flying 
qualities (see Military Standard 1797A sections: 4.5. 1.1 - Roll Mode, 4.5. 1.2 - Spiral 
Stability, 4.5. 1.3 - Coupled Roll-Spiral Oscillation, 4.6. 1.1 - Dynamic Lateral-Directional 
Response). This guidance is in the form of desired time constants, and frequency and 
damping specifications. 

Unfortunately, the Military Standard provides no direct guidance for choosing lateral- 
directional eigenvectors to yield desired flying qualities. Indirect guidance is provided in the 
form of lateral-directional modal coupling requirements. Two sections of Military Standard 
1797 A directly address lateral-directional coupling for relatively small amplitude rolling 
maneuvers (see Military Standard 1797A sections: 4.5. 1.4 - Roll Oscillations; and 4.6.2 - 
Yaw Axis Response to Roll Controller). In these sections, requirements are given placing 
limits on undesirable time responses due to control inputs. These requirements are based 
on time response parameters that can be measured in flight and were derived from flight 
data obtained from aircraft possessing conventional modal characteristics. The data base 
used to define the desired and adequate flying qualities boundaries (for high performance 
aircraft) is from Meeker and Hall, 1967. The data used to define this criteria is considered 
to be sparse. 

In addition to the Military Standard modal coupling criteria, some guidance is available 
from Costigan and Calico, 1989. The Costigan-Calico study correlated pilot handling 
qualities to the ratio of two elements of the Dutch roll eigenvector. Although this study did 
not lead to a design criteria, it does provide valuable pilot preference information for 
variations in the studied parameter. 
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The Meeker-Hall study, the Military Standard requirements for high performance 
aircraft performing precision tracking tasks, and the Costigan-Calico study are summarized 
in the following. 


Meeker-Hall Flight Test Study 

This study (Meeker and Hall 1967) examined the lateral-directional handling qualities 
for a selected range of dynamics for the fighter mission. Five groups of configurations with 
different combinations of lateral-directional dynamics were investigated. In three of the five 
groups, the effects of lateral-directional modal coupling were explored. In these groups, the 
roll, spiral and Dutch roll frequency and damping were held constant. Different modal 
couplings were achieved by varying the roll-to-sideslip ratio and bank angle-to-lateral stick 
transfer function numerator zeros. The roll-to-sideslip ratio l<J»/p \ dr is defined as the ratio of 
the amplitudes of the bank angle and sideslip time response envelopes of the Dutch roll 
mode, at any instant in time. Three roll-to-sideslip ratios were considered (1.5, 5-7, and 13). 
These groups (consisting of 58 configurations) were flight tested using a variable stability 
T-33 aircraft. One pilot evaluated each configuration using a series of tasks representative 
of the fighter mission (bank angle tracking, roll coordination, control in presence of 
disturbance, etc.). Subjective ratings were given in terms of Cooper ratings (Cooper 1957) 
and pilot comments. Modal characteristics, stability and control derivatives, transfer 
functions zeros, and pilot ratings for all three configurations are given in Tables 1-7. 

This study demonstrated that lateral-directional flying qualities are influenced by the 
relative location of the bank angle-to-lateral stick (or roll rate-to-lateral stick) transfer 
function numerator zeros with respect to the Dutch roll poles (equation (2.9) or (2.15)). 
The optimum pilot ratings occurred when the roll rate-to-lateral stick transfer function 
numerator zeros approximately canceled the Dutch roll poles. Configurations with zeros to 
the left of the Dutch roll pole were generally rated better than those with zeros to the right. 
In addition, the configurations with zeros in the lower left quadrant with respect to the Dutch 
roll pole showed less degradation in pilot rating as the zero was moved from its optimum 
location. (See Figure 1) 

For configurations with the lower roll-to-sideslip ratios the primary concern was the 
sideslip response that resulted from the lateral stick input rather than the roll response. 
These configurations have low coupling between the roll and sideslip responses and 
therefore the roll response is only slightly affected by large sideslip angles. For 
configurations with the medium roll-to-sideslip ratios the primary concern was the character 
of the roll response that resulted from the lateral stick input. Configurations with the larger 
roll-to-sideslip ratios (approximately 12 or greater along with a lightly damped Dutch roll 
mode) exhibited large rolling moments due to sideslip and were found to be unsatisfactory 
for the fighter mission. 

As a result of this study, specifying flying qualities criteria in terms of acceptable roll 
rate-to-lateral stick transfer function zero locations with respect to the Dutch roll pole was 
investigated This approach was found to have some major shortcomings. A primary 
shortcoming was the need to accurately determine the location of the zeros of the roll rate- 
to-lateral stick transfer function; this is difficult to measure. Industry preferred flying 
qualities criteria based on easily measured parameters. This lead to the development of the 
current time response parameter-based criteria in the Military Standard. 
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Military Standard Criteria 

As shown in the Meeker-Hall study, the existence of roll rate oscillations is directly 
traceable to the relative locations of the zeros and Dutch roll poles in the roll rate-to-lateral 
stick transfer function (equation (2.15)). When the complex roots cancel, the Dutch roll 
mode is not excited at all. When they do not cancel, there is coupling between the roll and 
sideslip responses. How this coupling is manifested depends upon the magnitude of the 
roll-to-sideslip ratio for the Dutch roll mode, l<t>/(3ld r • An approximation for the roll-to- 
sideslip ratio (Meeker and Hall 1967) is given by: 


( 


<t> 


b 

p 

dr 

Np 




1 + 


N p L r 
• 2 
b 


2 \ 


1 + ^ 


Nr. 


Vi 


(2.53) 


The Dutch roll contamination occurs primarily in yaw and sideslip if l<t>/fildr is low (less 
than approximately 1.5) or primarily in roll rate when l<|>/pldr is moderate-to-large (greater 
than 3.5 to 5). As l 4 »/( 3 ldr tends toward zero (Lp tends toward zero), the roll and sideslip 
responses become less coupled. 

In the Military Standard, pilot subjective flying qualities ratings are quantified in terms 
of Cooper-Harper ratings (Cooper and Harper 1969). The Cooper-Harper rating scale (and 
its predecessor the Cooper scale (Cooper 1957)), is a numerical scale from one to ten with 
one being the best rating and ten the worst (see Figure 2). In practice, Cooper-Harper 
ratings from one through three are referred to as "Level One", ratings from 4 through 6 as 
"Level Two", and seven through 9 as "Level Three". 


Roll Rate Oscillation Criteria 


The (pose I Pavg) parameter is directed at precision control of aircraft with moderate-to- 
high l<t>/pldr combined with marginally low Dutch roll damping. The ratio (p OS c ! Pavg) is a 
measure of the ratio of the oscillatory component of the roll rate to the average component 
of the roll rate following a step roll command (Chalk et al. 1969). This ratio is defined as 

Pose _ Pi + P3 ~ ' 2 l£1 (2 54) 

Pavg P\ + P3 + 2 P2 

for C,dr less than or equal to 0.2 and 


Pose _ P\ Pi 
Pavg Pi + Pi 


(2.55) 


for Q r greater than 0.2 where p\, P2, and p 3 are roll rates at the first, second, and third 
peaks; respectively. 

The values of (Pose / Pavg) that a pilot will accept are a function of the angular position 
of the zero relative to the Dutch roll pole in the roll rate-to-lateral stick transfer function. 
This angle will be referred to as 4 / i. Values of 'Ll for various zero locations are given in 
Figure 3. Because of the difficulty in directly measuring *Fi , the criteria is specified in 



14 


terms of the phase angle of the Dutch roll oscillation in sideslip (for a step input), *¥p (see 
equation (2.37)). This angle can be measured from the sideslip time response due to a step 
input. The angle 'Fi is directly related to the angle, *¥ p . For positive dihedral, this 
relationship is given by (Chalk et al. 1969) 


= ^1-270 (degrees) 


(2.56) 


This relationship is relatively independent of roll and spiral eigenvalue locations and 
holds for a wide range of stability derivatives. 

The (p osc / p avg ) criteria (for positive dihedral) is given in Figure 4. The Level One 
flying qualities boundary has a constant magnitude of 0.05 for 0 > 'Vp > -130 and -340 > 
'Vp > -360 degrees and a constant magnitude of 0.25 for -200 > *¥p > -270 degrees. The 
magnitude increases linearly from a magnitude of 0.05 at *¥p=- 130 degrees to a magnitude 
of 0.25 at *¥p =- 200 degrees. The magnitude decreases linearly from a magnitude of 0.25 
at 4^ =-270 degrees to a magnitude of 0.05 at *¥p =- 340 degrees. 

For all flying qualities levels, the change in bank angle must always be in the direction 
of the lateral stick control command. This requirement applies for step roll commands up to 
the magnitude which causes a 60 degree bank angle change in 1.7 7^ seconds, where Tj is 
the damped period of the Dutch roll eigenvalue. 





(2.57) 


The primary source of data upon which this (p osc / Pavg) requirement is based is the 
medium l<j)/pid r configurations of Meeker and Hall, 1967. 


Sideslip Excursion Criteria 

The (A/Jmax / kp ) parameter applies to sideslip excursions and is directed at aircraft with 
low-to-moderate l<|>/pidr • The term A/^ax is defined as the maximum sideslip excursion at 
the c.g. for a step roll command 


A/3 max =max(/J(0)-min(/J(0) forO<t<tp (2.58) 

where tp is equal to 2 seconds or one half period of the Dutch roll, whichever is greater. 
The term kp is defined as the ratio of “achieved roll performance” to “roll performance 
requirement” 


m 


<Preq 


t-t, 


req 


(2.59) 


where <p(t) is the bank angle at a specified period of time, t req and <p req is the bank angle 
requirement specified in the Military Standard (Chalk et al. 1969). For example, a <p req 
typically used for high performance aircraft is 60 degrees bank angle at one second. For 
this requirement, equation (2.59) reduces to 


kp = 


<P(t) 


60 


f=lsec 


(2.60) 
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where <|)(t) has units of degrees. 

Analysis of data revealed that the amount of sideslip that a pilot will tolerate is a function 
of the phase angle of the Dutch roll component of sideslip, *¥p (equation (2.56)). When 
the phase angle is such that /? is primarily adverse (out of the turn being rolled into), the 
pilot can tolerate a considerable amount of sideslip. When the phasing is such that p is 
primarily proverse (into the turn), the pilot can only tolerate a small amount of sideslip 
because of the difficulty of coordination. 

The (Afimax / kp ) requirement is given in Figure 5. As shown, the Level One flying 
qualities boundary has a constant magnitude of 2 for 0 > 'Vp > —130 and -340 > 'F/j > - 
360 degrees and a constant magnitude of 6 for -200 > ^p > —270 degrees. The magnitude 
increases linearly from a magnitude of 2 at *¥p =-130 degrees to a magnitude of 6 at p = 
—200 degrees. The magnitude decreases linearly from a magnitude of 6 at 'V p =-270 
degrees to a magnitude of 2 at ¥^=-340 degrees. 

This requirement applies for step roll control commands up to the magnitude that causes 
a 60 degree bank angle change within Td or 2 seconds, whichever is longer. The primary 
source of data from which the sideslip requirement for high performance aircraft evolved is 
the low i<|)/pidr (approximately 1.5) configurations of Meeker and Hall, 1967. In general the 
available data suggest that (A fimax / kp) is not as useful as (p osc / p avg ) when l<t>/pidr >3.5 to 
5.0. 


Costigan-Calico Flight Test Study 

The primary objective of this study (Costigan and Calico 1989) was to correlate pilot 
handling qualities to the magnitude of the roll-to-sideslip ratio for the Dutch roll mode, 
l<(>/pidr • Seven combinations of roll-to-sideslip ratio magnitude (l<|>/pidr ) (0,1.5, and, 3.0) 
and phase angle (Z(<f>/p)d r ) (0, 60, and 120 degrees) were tested (Table 8). These variations 
were achieved by varying the magnitude and phase of the ratio of the <J> and P elements of 
the Dutch roll eigenvector. System eigenvalues, and roll and spiral eigenvectors were set to 
desired values and not varied (Table 9). Control laws were designed using eigenspace 
assignment and flight tested on a YA-7D test aircraft. Three pilots evaluated each of these 
seven configurations using two Heads-Up-Display tracking tasks (yaw pointing and bank 
angle tracking) and an air-to-air task with a cooperative target. Pilot ratings were given in 
terms of Cooper-Harper ratings (Cooper and Harper 1969). Average Cooper-Harper 
ratings for the yaw pointing and bank angle tracking tasks are summarized in Table 10. 

Overall, the results showed little variation of the pilot ratings with l<t>/pid r and a 
preference for zero degree roll-to-sideslip phase angle over the larger phase angles. 
Costigan and Calico state that they believed the poor lateral stick dynamics of the YA-7D 
test aircraft degraded the lateral flying qualities ratings in all tasks and contributed to the 
small variations in pilot rating with k))/pidr • 


As shown, the Military Standard provides indirect guidance for choosing lateral- 
directional eigenvectors to yield desired flying qualities in the form of the roll rate 
oscillation and the sideslip excursion criteria. Using these criteria, it is not clear how to 
choose eigenvectors to achieve desired closed-loop flying qualities, or trade-off flying 
qualities for other important design requirements, such as achieving realizable gain 
magnitudes or desired system robustness. The next section addresses this shortcoming by 
presenting the development of guidelines for choosing lateral-directional eigenvectors to 
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yield desired flying qualities. This is done by developing relationships between the lateral- 
directional eigenvectors and the roll rate and sideslip transfer functions. Using these 
relationships, along with constraints imposed by system dynamics, key eigenvector elements 
are identified and guidelines for choosing values of these elements to yield desirable flying 
qualities developed. 


3.0 LATERAL-DIRECTIONAL EIGENVECTOR DESIGN GUIDELINES 

Eigenspace assignment methods allow designers to shape the closed-loop response by 
judicious choice of desired eigenvalues and eigenvectors. As shown earlier, the eigenvalues 
determine the time constant or the frequency and damping of each mode and the 
eigenvectors determine how much each mode of the system contributes to each output. 

When choosing desired closed-loop eigenvectors, the designer is faced with two 
challenges. First, the designer must choose which of the elements of the eigenvector matrix 
to specify. Using eigenspace assignment methods, for a system with n states, m controls, 
and / measurements, one has the freedom to place / eigenvalues to desired locations and m 
elements of their associated eigenvectors to desired values. Since for aircraft there are 
usually fewer controls than states, only a subset of the system eigenvectors can be exactly 
specified. Secondly, once the designer has chosen which elements to specify, he/she must 
decide what values to specify. Currently, no guidelines exist for choosing lateral-directional 
eigenvector elements to yield desirable flying qualities. Design guidelines would allow 
designers to perform trade-offs between flying qualities and other important design 
requirements, such as achieving realizable gain magnitudes or desired system robustness. 

This section addresses these two challenges by developing relationships between the 
system's eigenvectors and the roll rate and sideslip transfer functions. Using these 
relationships, along with constraints imposed by system dynamics, key eigenvector elements 
are identified and flying qualities guidelines for choosing appropriate values of these 
elements developed. 


Transfer Functions and Eigenvector Element Ratios 

Because eigenvectors can be scaled by an arbitrary constant, individual elements of an 
eigenvector are not unique. But, ratios of two elements of the same eigenvector are unique. 
Because of this, the eigenvector relationships and guidelines developed in this section will 
be stated in terms of these ratios. These ratios will be referred to as eigenvector element 
ratios. An eigenvector element ratio is equal to the ratio of the jc, and xj elements in the 
eigenvector associated with mode k and will be denoted by 

f \ 


V X j /mode k 

Eigenvector element ratios (also referred to as modal response ratios) can be expressed 
using any one of the n cofactors of the system's characteristic determinant (McRuer et al. 
1973). The eigenvector element ratio between two states x, and xj, and evaluated at mode k, 
is given by 
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V 


^ Xj 7mode k 




A gjjs) 

V A 4/W 


for q = 


(3.2) 


I s=X t 


where Aqi(s) and A qj(s) are the minors of the characteristic determinant A(j) of the system 
given by equation (2.1). 


Roll Rate Oscillations and Eigenvector Element Ratios 

The existence of roll rate oscillations is directly traceable to the relative locations of the 
zeros and Dutch roll poles in the p-to-8 st k transfer function (equation (2.15)). Equations 
for the frequency and real part of the Dutch roll pole are given by equations (2.4) and (2.5), 
respectively. Equations for the p-to-5 s tk zeros, as a function of the Dutch roll frequency 
and damping and lateral-directional coupling derivatives are determined as follows. 

Substituting equation (2.4) into equation (2.10) and making the assumptions of (2.12) 
yields 


<>i = (o 2 dr 


1- 


f K i \ 
Nt 


w 


N f>J) 


Subtracting equation (2.5) from equation (2.1 1) yields 


2(0^4, = 2(O d £dr 


f N. ^ 


K L 5J 


L r + 






N„- 




8_ 

Vo. 


(3.3) 


(3.4) 


As shown by these equations, the location of the zeros of the p-to-5 s tk transfer function with 
respect to the Dutch roll poles are primarily a function of the control coupling derivatives 
(N'§/ L's); and the stability coupling derivatives (L'p/N'p), (N' p - g/Vf), and L' r . 

The control derivatives are elements of the control matrix (B matrix of equation (2.22)) 
and are therefore independent of system eigenvectors. The ratio of control derivatives 
(N's/L's) can be adjusted by blending the lateral and directional control effectors (e.g. an 
aileron-rudder interconnect). The stability derivatives are elements of the state matrix (A 
matrix of equation (2.22)) and can be related to the system eigenvectors. This is done in the 
following. 


Solution for L' r 

The eigenvector element ratio (0//J) is given by applying equation (3.2) with i=2 and j= 1 


fi] 


f a 32 1 

-L r (s- Ypj + Lp 

UJ 


l A 3i ; 

s 2 — L n s — — L r 


(3.5) 


where q - 3 is chosen to yield a desired solution form. Evaluating this ratio at 5 = X dr , and 

2 

recognizing that both YpL' r « L'p and ( g/Vo)L' r «(O dr yields 
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f </> ^ _ ~ L r^dr + L B 

R I “ ; 2 _ t' x 

V P J dr A'dr ^p^dr 

The term (<t>/p)d r is the ratio of the bank angle and sideslip elements in the Dutch roll 
eigenvector. The phase angle of equation (3.6) is given by 


Z ^ = -atan -4 


L r (O d ri I"! 


+ Cdr^dr^r 


+ atan 


®dr\^ £dr (^7 ? + ^Cdr^dr ) 
^dri^Cdr ~ l) Cdr^dr^p 


= atan 


^ dr yl Cdr 2 Cdr®dr^f) ^dr ^-r j 

^dr{^Cdr ~ Cdr^dr^r Cdr^drLpLp + G>d r L r Lp 


Solving this equation for L' r yields 


L r = - 


( L p + ^dr^dr )<X>dr ^~Cdr “ (®dr( 2 Cdr ~ l ) + Cdr a dr^p) tan Z ~ 


^dr Cdr \^P ^dr^dr j ^an •*— „ 

VP 


This phase angle Z(4>/p)d r is related to the phase angle Z(p/f$)dr by the following relation 


The phase angle Z(p/p)dr is a discriminator of positive and negative dihedral (Chalk et al. 
1969). Positive dihedral corresponds to 45° < Z(p/($)dr < 225°. For a stable complex 
Dutch roll pole 


Z X dr = 180° - <2cos(£^ r ) 

Therefore, for positive dihedral (and a stable complex Dutch roll pole) 


-135° + acos(Cdr) < Z^J < 45° + acos(C dr ) 


(3.10) 


(3.11) 


where 0° < acos(& r ) < 90°. Results from Costigan and Calico (1989) show a pilot 
preference for Z(<})/P)dr =0. Choosing Z(4>/P)<j r =0 provides positive dihedral for any stable 
complex Dutch roll pole and simplifies the solution for L' r . For Z(<t>/P)d r =0, equation 
(3.8) reduces to 
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A relationship between X ro u , L' p , and coupling derivatives is given by equation (2.6). 
Solving this for L' p yields 


l p = x ro u + -4 \n'-£- 

p rotl No I p 


Substituting this equation for L' p (equation (3.13)) into (3.12) yields 


L r = ~ ~r Koll + TT N P~T7 + 2 Cdr<°dr 
No Na Vn J 


(3.13) 


(3.14) 


Solution for (N’ p - g/Vo) 

The eigenvector element ratio (fi/p) is given by applying equation (3.2) 


g) (A3,) 


(3.15) 


\P) L p -L r (s-Y p ) 

Evaluating this ratio at s = X ro u and recognizing that both YpL' r «L'p and 
(g/Vo)L'r/X ro u « L'p yields 

'£} s (3.16) 

\P ) roll Lp~Koll L r 

The term (p/p) r oll is the ratio of the sideslip and roll rate elements in the roll mode 
eigenvector. Substituting equation (3. 1 3) for L' p into (3.16) yields 


iU'-i 




£ _ 

p J roll L p~Koll L \ 


(3.17) 


2 1 

Substituting equation (3.14) for L' r into (3.17) and applying (Oj r = Np yields 


£ 

p Kou 


- N p --2- 

{ p v 0 


03 dr + Koll Koll + Kdr^dr + ~7~ II N p~~ 


(3.18) 
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Solving for ( N' p -g/Vo ) yields 



(3.19) 


Solution for (L'p/N'p) 

Substituting equation (3. 14) for L' r into (3.6) and applying 6)j r = Np yields 


\P'dr 


J >1 

b 

[Lp + 2^ dr ^dr)^dr + ®dr 


1 

4 1 

1 

i 


(3.20) 


By noting that the term in square brackets is equal to - 1 (and that for Z(< t>/P)<i r =0 the ratio 
(4>/P)dr is a real number) equation (3.20) reduces to (for positive dihedral) 




L P 


A r, 


P. 


(3.21) 


By substituting equations (3.14), (3.19), and (3.21) into (3.3) and (3.4); one can now obtain 
expressions for the frequency and damping of the transfer function zeros as a function of 
0)dr, Cdr , &roli . (N'§/ L's), and eigenvector element ratios. This is done in the following 
steps. 


Substituting equation (3.21) into (3.3) yields 


<4 = <°dr 


1 + 





(3.22) 


Substituting equation (3. 14) for L' r into (3.4) yields 


— 2 COdr^dr 




Ns 


N, 


Koii + 


L P 




PA 1 * 


N, 


P 


f ■ g \ 

N - — 

P 




Vr 


0J 


+ 


f f \ 

b 

\ N p j 


( 


N, 



+ -^dr^dr 

( 3 . 23 ) 



Rearranging terms yields 


2 ^<f)C<p — 2 ®dr^dr 


( j \ 

b 

K n 'p 


7^ ( Koll + ^dr^dr) ~ [ N p ~ “1 

\ L SJ v v o) 


f r f' V »r \ 


1- 




N, 


P) 


N* 

vbj 


M 


J J 


(3.24) 


Now substituting equation (3.19) for ( N' p -g/Vo ) into (3.24) yields 
— 2 ®drCdr ~ 


' b_ 

KJ 


{Koll +2 Cdr (0 dr) + {— 

i L s) yp 


P ^ {^roll + ^bollCdr^dr + ^dr) 


'roll 


1 + A, 


ro// 




f ,■ Y\ 


yPJroll{ N P 


± 

Nr. 


1- 


V^A 4 JJ 


J) 


(3.25) 


Finally, substituting equation (3.21) for (L'p/N’p) into equation (3.25) yields the desired 
relationship 


2 ®<f>C<p — 2 ®drCdr 4 " 


fi) 

^Pjdr 


f u' A 

ybj 


broil + ZCdr^dr) 


fi \ [fi-oll + 2 Koll Cdr^dr + 03 dr )( f <p 

f?l ) 1 u 


dr 


Nr 




\ L 5 )) 


PJrollVPjdr, 


(3.26) 


Equations (3.22) and (3.26) provide equations for (Q<p and as a function of system 
eigenvalues, eigenvector element ratios ( l<|>/pidr and (p/P) ro ii ) and the adverse yaw ratio 
(N's /L's). The relationship between (p/P) r oll and (N'$ /L'p) and co# and ^ is 
demonstrated in Figure 6 for l<t>/pldr = 5 and eigenvalues X ro u = -2.5, A^ r = ( (0d r = 2.0, Q r 
= 0.1). In this figure, lines of constant (D<p are given by solid lines (1.5, 2, 2.5, 3 (rad/sec)) 
and lines of constant are given by dashed lines (0, 0.4, 0.7, 1). 
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As can be seen from equation (3.26), the value of the ratio l<(>/pidr has a strong effect on 
the transfer function zero locations. For l<t>/pidr = 0, the locations of the transfer function 
zeros are not affected by variations in (p/(3) ro li and (N'g/L's). As l<t>/pl dr increases, the 
zeros become more sensitive to variations in (p/P) ro ll and (Nf/Ls ). 


Sideslip Response and Eigenvector Element Ratios 

The sideslip-to-lateral stick transfer function can be written as (equation (2.19)) 


-N* 


P 


s + —L n + - 


Ns 


N p - — 
P V f 


Vi 


0 ))) 


^stk (5 ^roll ^Cdr®dr s + 

Rearranging terms in the numerator yields 



^stk (s boll'){ s ~br^dr s j 


(3.27) 


(3.28) 


By substituting equations (3.13), (3.19), and (3.21) into (3.28); one can now obtain an 
expression for the sideslip-to-lateral stick transfer function zero as a function of system 
eigenvalues, eigenvector element ratios, and adverse yaw ratio. This is done in the following 
steps. 


Substituting equation (3. 13) for L' p into equation (3.28) yields 


—Ls 


jL 

%tk 


N. 


■r~(s-X r - 0 ii)+ N p ~~ 


i-Ofi. 


( / ^Y\ 


IP 


N 


ej)) 


( 5 “ Kolltf 2 + 2 br<°dr s + 


(3.29) 


Now substituting equation (3.19) for (N' p -g/Vo) into (3.29) yields 
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(3.30) 
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Finally, substituting equation (3.21) for (L'p / N'p) into equation (3.30) yields the desired 
relationship 



—Lz 


Koii )- 1 


i+- 
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] i^roll + ^A-rol&dr^dr + ®dr ) 

yPJroll 


L S \PJdr) 


1 -Koll - 


{ s ^ 'roll )(■* ®drj 


' roll 


^P'drJ 


(3.31) 


Equation (3.31) provides an equation for the sideslip response as a function of system 
eigenvalues, eigenvector element ratios ( l(])/pidr and (p/fi)roll ) and the adverse yaw ratio 
(N's/L’s). 


Additional Eigenvector Element Ratio Relationships 

This section presents additional eigenvector element ratio relationships imposed by the 
system dynamics. 


Roll Rate-to-Bank Angle Ratios 

The roll rate-to-bank angle ratios for each mode can be determined from the kinematic 
relationship between roll rate and bank angle 


_P_ = J±_ 

&stk $stk 


Evaluating this equation at each mode yields the following relationships 


f =V; 

‘fJdr 
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(3.32) 


(3.33) 


Yaw Rate -to-Side slip Ratio in the Dutch Roll Mode 

The yaw-to-sideslip ratio in the Dutch roll mode can be determined from the sideforce 
equation of (2.1). This is given by (making the assumptions of (2.12)) 


L_A- 
P V 0 


vPj 


+ (s-Yp) = 0 


(3.34) 


Evaluating this equation at s = Xdr > and solving for (r//J)dr yields 
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(3.35) 
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In general, both Yp and (g/Vo)(ty P)dr are very small compared to ^j r , therefore 


(3.36) 


Yaw Rate-to-Roll Rate in Roll Mode 

The yaw-to-roll rate ratio in the roll mode can also be determined from the sideforce 
equation of (3.34). Evaluating this equation at s=X ro u , and solving for {r/fi) ro u yields 


4] =-(A ro „-y fl ) + -«-W 

PJnX ' P ’ 
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'r' 


Multiplying both sides by (p/p) r oll and applying equation (3.33) yields 

yPKoll yPfrM 1 

In general, Yp is very small compared to , therefore 
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(3.38) 


(3.39) 


Yaw Rate-to-Bank Angle Ratio in the Spiral Mode 


The yaw-to-bank angle ratio in the spiral mode can be determined from the sideforce 
equation of (3.34). Evaluating this equation at s=X spr i , and solving for ( r/p) spr i yields 




ypjspn 


- -Kprl + Y p + -~ 


t 

( 3 , 
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(3.40) 


Multiplying both sides by (p/<f>) spr i and assuming (p/<t>) spr [ (Yp - X sprl ) is small compared 
to ig/Vo) yields 
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_ 8 


V 


tJsprl V 0 


(3.41) 


Specifying the Desired Eigenvectors in terms of Eigenvector Element Ratios 

Desired eigenvector elements will be chosen by considering both the relationships 
between the eigenvector element ratios and transfer functions, and constraints imposed by 
system dynamics. Desired eigenvector elements will be specified for the lateral-directional 
state vector given by 


x = { fipr<f> }T 


(3-42) 
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For lateral-directional design, usually only two controls are available - roll moment and 
yaw moment (the direct sideforce generated by conventional controls is usually small). This 
development will assume that four measurements (or full-state feedback) are available. 
Therefore, using eigenspace assignment methods one can exactly place all four lateral- 
directional eigenvalues and only two elements of each associated eigenvector. 

Of course, more than two elements in each desired eigenvector can be specified. This 
results in closed-loop eigenvectors that are as close as possible in a least squares sense to 
the desired eigenvectors. In this authors experience, for fixed desired eigenvalues, the 
resulting closed-loop eigenvectors are usually a poor fit to the desired eigenvectors. 
(Although not considered in this work, eigenvector fit can be improved by allowing 
variations in the desired eigenvalues within a certain region). Therefore, this work will only 
consider specifying two elements of each eigenvector. 

Because eigenvectors can be scaled by an arbitrary constant, one element of each 
eigenvector will be specified to be unity to ensure that the eigenvectors are unique. In the 
following desired eigenvectors, a * indicates that the element is not specified and therefore 
not weighted in the cost function. 


Spiral Eigenvector 

The spiral mode is a first order mode with a long time constant. Classically, the spiral 
mode is dominant in bank angle and almost nonexistent in sideslip. Therefore, the bank 
angle element will be chosen to be unity. The eigenvector element ratios available to be 
specified for this eigenvector are (p/<p ), and (r/0). Having a very small spiral mode 

contribution to sideslip is desirable because it results in coordinated banking and turning. 
This can be achieved by choosing the to be zero. This results in Rfrsprl = 0- The two 
remaining ratios (r/cj>) and (p/<p) are constrained by the system dynamics. The ratio of the p 
and 0 elements is constrained to be equal to the spiral eigenvalue (equation (3.33)) and the 
ratio of the r and <f> elements is constrained to be approximately equal to (g/Vo) (equation 
(3.41)). Therefore, the desired spiral mode eigenvector is specified to be 

«V( = 1° * * ‘f (3-43) 


Roll Eigenvector 

The roll mode is a first order mode with a relatively short time constant. Classically, the 
roll mode is dominant in roll rate therefore this element will be chosen to be unity. 
Therefore, the eigenvector element ratios available to be specified for this eigenvector are 
(/fy>), (r/p), and (<p/p). As was shown by equation (3.25), the ratio of the /? and p elements, 
(fi/p)roll - effects the cancellation of the Dutch roll mode in the roll response. The two 
remaining ratios, (r/p) and (<p/p), are constrained by the system dynamics. The ratio of the <p 
and p elements is equal to the inverse of the roll eigenvalue (equation (3.33)) and the ratio of 
the r and p elements is a function of the roll eigenvalue and (fi/p) (equation (3.39)). 
Therefore, the roll mode eigenvector is specified to be 
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Jroll 


—roll = 


T 


(3.44) 
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Dutch Roll Eigenvector 

The Dutch roll mode is a lightly damped oscillatory mode. For this eigenvector, the 
sideslip element will be chosen to be unit magnitude with zero phase. Therefore, the 
eigenvector element ratios available to be specified for this eigenvector are (p/ft), (r/ft), and 
Specifying the ratio of the ft and 0 elements in the Dutch roll eigenvector determines 
the l<t>/|3ldr ratio. Results from Costigan and Calico (1989) show a pilot preference for Z 
(<t>/P)dr = 0. The two remaining ratios, (p/ft) and (r/ft) are constrained by the system 
dynamics. The ratio of the r and ft elements is constrained by the system dynamics to be 
approximately equal to the negative of the Dutch roll eigenvalue (equation (3.36)). The ratio 
(p/ft) can be written as the product of the (<///?) and (p/<j>) ratios. Because the ratio of the p 
and 0 elements is constrained to be equal to the Dutch roll eigenvalue (equation (3.33)), 
once the (<yft) ratio is specified then the (p/ft) ratio is also specified. Therefore, the Dutch 
roll mode eigenvector is specified to be 
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<p 
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where ( • , • ) denotes (magnitude, phase (degrees)). 


Eigenvector Element Ratio Design Guidelines 

In the previous section, desired lateral-directional eigenvectors were specified in terms of 
eigenvector element ratios ( l<t>/|3ldr and (p/ft) ro ll ), and the adverse yaw ratio (N's / L\ 5 ). In 
this section, relationships between the Military Standard lateral-directional coupling criteria 
and these parameters will be developed. Using these relationships, the Military Standard 
Level One flying qualities boundaries can be translated into guidelines on eigenvector 
element ratios and the adverse yaw ratio. 


Roll Rate Oscillation Criteria and Eigenvector Element Ratios 

The Dutch roll contamination occurs primarily in roll rate when l<|>/{3ldr is moderate-to- 
large. For these configurations, the Dutch roll contamination can be quantified in the time 
domain by the ratio (p osc / p avg ) (equations (2.54-55)). This ratio can be expressed as a 
function of l<j)/pidr , (p/ ft) roll > and (N's/ L's ) by defining it in terms of the system residues. 
This is done in the following. 

The roll rate due to a unit step input in lateral stick expressed in terms of system 
residues and eigenvalues is given by (because there is only one input, the third subscript on 
the R 's has been omitted) 
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where (Churchhill et al. 1976) 
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x = L dKr(Kr Z<p){^dr ty) 

P (irfr ~ Kprl )(Kr ~ Koll ){^dr ~ A dr ) 

and z<p = -(Oq^q + and x denotes the complex conjugate of x. The p-to-S st k 

transfer function numerators z<p can be expressed in terms of system eigenvalues, l<t>/pldr , 
(p/P)roll , and (N's/L's) as 
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where x denotes an approximation to x. An approximation to p(t) expressed in terms of Zq 
is given by 
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The ratio (p osc /p avg ) can be expressed as a function of p(t) by 

Pose ~ Pose _ Pi + Pi ~ 2 P2 
Pavg Pavg Pl + P3+ 2p2 


for Cdr less than or equal to 0.2 and 


Pose ~ Pose _ P\ P2 
Pavg Pavg Pi + P2 


(3.55) 


(3.56) 


(3.57) 


(3.58) 


f° r Cdr greater than 0.2 where p x , p 2 , and p 3 are values of p(t) at the first, second, and 
third peaks; respectively. 


Equations (3.57-58) (along with equations (3.50-56)) provide expressions for 
(Pose/ Pavg ) as a function of system eigenvalues, 10/pldr , (p/fyroll , and (N's / L’$). Using 
these equations (p osc f p avg ) can be calculated in the following way: 

1) Choose system eigenvalues and values of l<))/pidr , (p/P)roll , and ( N's/L's ). Note 
that this development requires (1 + (N's/L's) l<t>/pldr )>0. 

2) Calculate 0 )q and ^ using equations (3.51-52) and form 1$. 

3) Calculate R pspr i, R pjo lb and Rp,dr using equations (3.54-56) and form p(t) . 

4) Generate step time response using p(t) (equation (3.53)). 

5) Pick off peaks from pit) step time response (p\, p 2 , and p^). 

6) Calculate ( Pose /Pavg ) using equation (3.57) or (3.58). 


Sideslip Excursion Criteria and Eigenvector Element Ratios 

The Dutch roll contamination occurs primarily in sideslip if l<j>/pi dr is low. The Dutch 
roll contamination can be quantified in the time domain by the ratio (A/l^ax / kp ) (equations 
(2.58-59)). This parameter can be expressed as a function of !<}>/pidr , (p/ ft) roll , and 
C N's/L's ) by defining it in terms of the system residues. This is done in the following. 

The sideslip due to a step input in lateral stick expressed in terms of system residues 
and eigenvalues is given by (because there is only one input, the third subscript on the R 's 
has been omitted) 
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Making the assumptions of (2.12) and by specifying the desired spiral eigenvector as 
defined in (3.43), equation (3.59) reduces to 
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An approximation to 0(r) expressed in terms of l<|)/pidr , (p/fi)roll , and (N'§/L's) is given by 
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and 5: denotes an approximation to x. The ratio (A/Jmax / kp ) can be expressed as a 
function of pit) by 
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where 
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where Rp S p r i , R pro n, and are defined by equations (3.54-56). 



Equations (3.69-71) (along with equations (3.65-3.68) and (3.72-74)) provide 
expressions for (A/) max / kp) as a function of system eigenvalues, eigenvector element 

ratios ( l<J>/pldr and (p/fi)roll )> an d (N'$/ L'§) . Using these equations (A/) max /kp) can be 
calculated in the following way: 

1) Choose system eigenvalues and values of l<|>/pldr . (p/p)roll > and (N'$ / L'§) Note 
that this development requires (1 + (N'§/ L'$) l(J)/pid r )>0. 

2) Calculate /3q , Rp ro u , and Rp^ r (using equations (3.66-68)) and form P(t). 

3) Generate step time response using j3(t ) (equation (3.65)). 

4) Calculate A/) max (equation (3.70)) from (3(t) step time response. 

5) Calculate 0 q, R^^ spr i, ra //, and R$j r (using equations (3.73-74)) and form 0(t). 

6) Calculate kp (equation (3.71)) from (equation (3.72)). 

7) Calculate (A/J max /kp) using equation (3.69). 


Phase Angle of Dutch Roll Component of Sideslip and Eigenvector Element Ratios 

The values of (p 0 sc / p a vg ) and (A/Jmax ) that a pilot will accept are a function of the 
phase angle of the Dutch roll component of sideslip, 'Vp . This angle is directly related to 
the angular position of the zero relative to the Dutch roll pole in the roll rate-to-lateral stick 
transfer function 4 / i. This angle is given by 


f 

'Pi = atan 




^dr ^dr ®^^1 4 ^ 

~Cdr®dr + C<j>^<t> j 


(3.75) 


As shown earlier (equation (2.56)), for positive dihedral, 'Vp and *Fi can be related by 
(Chalk et al. 1969) 


Vp a - 270 (degrees) 


(3.76) 


An approximation to 'Vp expressed in terms of system eigenvalues, l<J>/pidr » (p/fi)roll > and 
(N's/L's) is given by 
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where and ^ are defined by equations (3.51) and (3.52). 
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Equation (3.77) provides an expression for ^ as a function of system eigenvalues, 

l0/pidr > (p/p)roll » and (N's/ L's). Using these equations can be calculated in the 
following way: 

1 ) Choose system eigenvalues and values of t4>/pidr . (p/P)roll , and (N's /L's). Note 
that this development requires (1 + (N's/ L's) l<J>/pid r )>0 

2) Calculate d ) and ^ using equations (3.51-52). 

3) Calculate 4^ using equation (3.77). 

This development assumes values of l<t>/(3l dr , (p/p) ro ll , and (N's/ L's) have been chosen 
that yield complex conjugate transfer function zeros (i.e. -lc^cl and d)^ > 0 ). 


Generating Eigenspace Flying Qualities Guidelines 

Eigenspace flying qualities guidelines for choosing (p/(5) r oll and (N's/ L's) for a given 
value of l<j)/|3ldr are determined in the following way: 

1) Choose value of l<t>/pi d r and desired eigenvalues (roll and spiral eigenvalues at Level 
One locations). 

2 ) Evaluate ( Posc/Pavg )> (A/? max Ikp), and 4^ over desired range of (p/(5) ro ll and 
N 5 /L 5 . 

3) The eigenspace roll rate oscillation guideline is based on the Military Standard roll 
rate oscillation criteria. This guideline can be determined by overlaying the ( Posc/Pavg ) 

and 4 ^ data and translating the Level One roll rate oscillation boundary into boundaries on 
(p/P)roll and (N's/ L's) . 

4) The eigenspace sideslip excursion guideline is based on the Military Standard 
sideslip excursion criteria. This guideline can be determined by overlaying the 

(A^ma x/kp) and 4^ data and translating the Level One sideslip excursion boundary into 
boundaries on (p/f}) r0 u and (N’s/ L's ) . 

The appropriate guideline to use is a function of the value of l<J>/(3ldr • For low |(J)/p| dr 
(less than approximately 1.5), the eigenspace sideslip excursion guideline should be used. 
For moderate-to-high l<t>/pi dr (greater than approximately 5), the eigenspace roll rate 
oscillation guideline should be used. For intermediate values of l<J>/pid r , both guidelines 
must be satisfied to meet Level One flying qualities. A composite eigenspace roll rate 
oscillation/sideslip excursion guideline can be obtained by overlaying the guidelines from 
these two criteria. 


Eigenspace Flying Qualities Guideline Examples 

Example guidelines are generated for a moderate-to-large, intermediate, and a low value 
of l4>/pl dr • These guidelines are developed for values of l<|>/pidr , (p/fi)roll , and (N’s /L's) 

that yield complex conjugate transfer function zeros in the left half plane (i.e. Ck^cl and 

0). For all of these cases the eigenvalues are set to the following values: 
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X spr i = -0.005, X roii = -2.5, Xdr- ( C0dr = 2.0 (rad/sec), Q r = 0.1). The values of l<f>/pidr 
considered are: kjj/pidr = 10, l<f>/pidr = 5, and l<|>/pldr = 1. 

Example One - \<t>/j3\dr —10 

The roll rate oscillation criteria is the suggested primary criteria for moderate-to-large 
l<t>/Pldr • The first step in generating the eigenspace roll rate oscillation guideline is to 

evaluate (p osc /p avg ) 311(1 over the desired range of (p/ ft) roll and ( N's /L's ). Figure 7 
shows how lines of constant ( Posc/Pavg ) are a function of (p/(3) ro ll and (AT 5 / L's) for 
l<f>/Pldr = 10. This figure shows contour lines for ( Posc/Pavg ) equal to 0.05 and 0.25. The 
magnitude of ( Posc/Pavg ) 1S directly related to the cancellation of the Dutch roll pole in the 
roll rate-to-lateral stick transfer function. At (p/(3) ro ii =0 and (Nf/L$ )=0, the Dutch roll 
pole is canceled and ( Posc/Pavg ) equals zero. As (p/{3) ro ll and (N's / L's) increase in 
magnitude, there is an increase in the Dutch roll modes contribution to the roll rate response 
and therefore ( Posc/Pavg ) increases in magnitude. 

Figure 8 shows how lines of constant 4^ are a function of (p/(3) r oli and (N's/ L's) for 
a value of I4>/Pldr = 10. This figure shows contour lines for 4^j equal to 0, -130, —200, 

-270, and -340 degrees. Since 4yj is directly related to the angular position of the zero 
relative to the Dutch roll pole in the roll rate-to-lateral stick transfer function, the constant 
‘F/j lines radiate from the point (p/f3) ro ll =0, (N's/ L's)=0. By overlaying the ( Posc/Pavg ) 

and 4^ contour data, the Military Standard Level One roll rate oscillation criteria can be 
translated into boundaries on (p/P) ro ll and (N's / L's). This is done in Figure 9. This 
figure shows how the four line segments that compose the Military Standard criteria map 
into the eigenspace guideline. The result of this mapping is given by the solid black line. 
Values of (p/ (3) roll and (N's / L's) that lie inside this boundary meet this criteria. This 

guideline is presented again in Figure 10 without the ( P osc /Pavg ) and ’P/j contours. 
Example Two — \<p/(3\dr = 5 

For intermediate values of l<t>/pidr , both guidelines must be satisfied to meet Level One 
flying qualities. The eigenspace guidelines are generated by first evaluating ( Posc/Pavg )» 

( ^Anax Jkp ) > and 4^ over the desired range of (p/P) ro ii and (N's /L's). The eigenspace 

roll rate oscillation guideline is obtained by overlaying the ( P OS c/Pavg ) and 4 # « data and 

translating the Military Standard Level One roll rate oscillation boundary into boundaries on 
(p/ (3) roll and (N's/ L's ) . This guideline is given by the solid black line in Figure 11. The 

eigenspace sideslip excursion guideline is generated by overlaying the (A/3 max /£q) and 

4^ data and translating the Military Standard Level One sideslip excursion boundary into 

boundaries on (p/f3) ro ll and (N's /L's). This guideline is given by the solid gray line in 
Figure 12. 

Both of these criteria must be satisfied to meet Level One flying qualities. The 
composite eigenspace guideline is obtained by overlaying these two criteria. This is done in 
Figure 13. The eigenspace roll rate oscillation guideline is given by the solid black line and 
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the eigenspace sideslip excursion guideline is given by the solid gray line. Values of 
(p/fi)roll and (N's/L's) that lie inside both of these boundaries will yield Level One flying 
qualities. As can be seen, for this value of l<J)/(3ldr both guidelines are approximately the 
same size, but over most of the boundary the roll rate oscillation guideline is the more 
restrictive criteria. 

As can be seen by comparing Figures 10 and 11, the eigenspace roll rate oscillation 
guideline is approximately the same shape as for l<f>/pidr =10 but considerably larger in size. 
This demonstrates the sensitivity of this guideline to the value of l<|)/pldr • As l<t>/J3ldr 
decreases, the system's flying qualities become less sensitive to variations in (p//5) ro ll and 
(N's/L's). 

Example Three — \(p/p\dr — 1 

The sideslip excursion criteria is the suggested primary criteria for low l<()/pidr- The 
eigenspace sideslip excursion guideline is generated by first evaluating (A(3 max j kp) and 
Pp over the desired range of (p/ ft) roll and (N's / L's). The guideline is obtained by 

overlaying the (A/3 max /kp) and Pp data and translating the Military Standard Level One 
sideslip excursion boundary into boundaries on (p/(3) ro ll and (N's/ L’s). This guideline is 
given by the solid gray line in Figure 14. Values of (p/f}) ro ll and ( N's /L's ) that lie inside 
this boundary meet this criteria. 

As can be seen by comparing Figures 12 and 14, this guideline is approximately the 
same size as for l<|>/pidr = 5. This guideline is considerably less sensitive to variations in the 
value of l<t>/pidr than the roll rate oscillation guideline. 


4.0 CONCLUDING REMARKS 

This report presents the development of lateral-directional flying qualities guidelines 
with application to eigenspace assignment methods. These guidelines will assist designers 
in choosing eigenvectors to achieve desired closed-loop flying qualities or performing trade- 
offs between flying qualities and other important design requirements, such as, achieving 
realizable gain magnitudes or desired system robustness. This has been accomplished by 
developing relationships between the system's eigenvectors and the roll rate and sideslip 
transfer functions. Using these relationships, along with constraints imposed by system 
dynamics, key eigenvector elements have been identified and guidelines for choosing values 
of these elements to yield desirable flying qualities developed. 

Two guidelines are developed - one for low l<|>/pidr and one for moderate-to-high 
l^/pidr • These flying qualities guidelines are based upon the Military Standard lateral- 
directional coupling criteria for high performance aircraft. The low l4»/|3ldr eigenspace 
guideline is based on the sideslip excursion criteria. The high l<j>/pi<j r eigenspace guideline 
is based on the roll rate oscillation criteria. For intermediate values of l<)>/|3ldr , both 
guidelines must be satisfied to meet desired flying qualities. A composite eigenspace roll 
rate oscillation/sideslip excursion guideline is obtained by overlaying the guidelines from 
these two criteria. 

Example guidelines are generated for a large, an intermediate, and low value of l<f>/pidr- 
For all of these cases the eigenvalues are set to fixed values and are not varied. In these 
examples it was shown that the value of the ratio l<t>/pid r has a strong effect on the 



eigenspace roll rate oscillation guideline, whereas the eigenspace sideslip excursion 
guideline is relatively insensitive to variations in l<t>/pidr- 

Piloted simulation flying qualities experiments are planned to validate and refine these 
guidelines. 
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Appendix - Primary Lateral-Directional Coupling Derivatives 

A coupling derivative is one in which a motion, angle, or control input about one axis 
imparts a moment about an orthogonal axis. The primary lateral-directional coupling 
derivatives are: roll moment due to sideslip angle Lp , roll moment due to yaw rate L r , yaw 
moment due to roll rate N p , and yaw moment due to lateral controls N§. 


Adverse Yaw due to Ailerons Ns 

In an aircraft in which ailerons are the primary lateral control effector, a right roll control 
input results in the left aileron down and right aileron up. This leads to a rolling moment 
due to more lift on the left wing and less on the right. More lift on the left wing results in 
more drag on the left wing. Therefore, there is also a yaw moment applied to the aircraft. 
For ailerons this yaw moment is opposite to the turn (adverse). 


Dihedral Effect Lp 

Dihedral of the open-loop airframe is primarily affected by wing location (mid-wing, 
high-wing, etc.), dihedral angle, and sweep. When an aircraft with positive dihedral 
encounters a positive sideslip it will tend to roll to the left because the right wing will see a 
higher angle-of-attack. An aircraft with positive dihedral effect will cause an aircraft to roll 
away from the sideslip. A negative value of Lp is positive dihedral. 


Roll Moment due to Yaw Rate L r 

When an aircraft rotates to the right, the left wing will see an increase in forward velocity 
(and the right a decrease) due to the rotation. This velocity change results in lift changes 
that cause a roll moment in the direction of the yaw rate. In addition, the yaw rate generates 
a lateral velocity change at the tail. This results in a side force at the tail (which is usually 
above the roll axis) causing a roll moment. 


Yaw Moment due to Roll Rate N p 

The primary contribution to N p is due to the wing. This derivative is a component of 
adverse yaw. When an aircraft rolls to the left, an angle-of-attack increment is generated 
due to the roll rate. This angle-of-attack increment increases the lift on the downward wing 
(and decreases on the upward wing) and results in the lift vector being tilted forward. This 
results in a yaw moment usually opposite to the roll. 


Primed Derivatives 

The primed derivatives are defined by 

j _ + . xj _ Nj +(l xz /I z )L; 

■ -(4Aa) ' ' i -(i 2 J'j z ) 

where the subscript i denotes a motion or input quantity. 


(A.l) 
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Figure 1 - Areas of acceptable zero locations p-to-lateral stick transfer function. 
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Figure 2 - Cooper-Harper handling qualities rating scale (Cooper and Harper 1969). 
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Figure 4 - (p osc / p a vg) Roll Oscillation Criteria for Positive Dihedral. 
(From Military Standard 1797 A). 
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Figure 5 - (A fimax / kp) Sideslip Excursion Criteria for Positive Dihedral. 
(From Military Standard 1797 A). 
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(N't/L's) 



Figure 1 1 - Level One Eigenspace Roll Rate Oscillation Guideline for l^/pidr = 5. 
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(N'z/L'g) 



Figure 13 - Overlay of Eigenspace Roll Rate Oscillation (black line) and Eigenspace 
Sideslip Excursion (gray line) Flying Qualities Guidelines for l<|)/pidr = 5. 
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Figure 14 - Level One Eigenspace Sideslip Excursion Guideline for l<t>/(3ldr = 1- 
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AB-1 

AB-2 

AB-3 

( 0 d r (rad/sec) 

2.36 

2.51 

2.55 

Cdr 

0.12 

0.10 

0.11 

WPldr 

1.43 

1.59 

1.34 

Z((j)/p)dr(deg) 

49.0 

44.0 

47.0 

Troll (sec) 

0.4 

0.4 

0.4 

Tspiral (sec) 

-91.0 

995 

OO 


Table 1 - Modal Data for Meeker-Hall Group AB. 
(Data from Meeker and Hall 1967) 


Sub- 

Group 


AB-1 


AB-2 


AB-3 


Config. 

City/COdr 

€00 

(rad/sec) 

Co 

L'S 

(deg/sec 2 )/in 

N's/L’s 

N 'p 

(sec) 

Pilot 

Rating 

a 

0.85 

2.01 

0.05 

66 

-0.166 

0.211 

7 

b 

0.88 

2.07 

0.06 

71 

-0.142 

X 

X 

7 

c 

0.91 

2.15 

0.06 

73 

- 0.111 

7 

d 

0.96 

2.27 

0.07 

72 

-0.066 

X 

6.5 

e 

1.00 

2.36 

0.08 

91 

-0.027 

X 

7 

f 

1.03 

2.42 

0.08 

70 

0.0 

0.211 

8 

a 

0.78 

1.97 

0.09 

75 

-0.159 

-0.0125 

7 

b 

0.85 

2.12 

0.10 

127 

- 0.111 

X 

4 

c 

0.93 

2.34 

0.10 

153 

-0.038 

X 

3 

d 

0.96 

2.40 

0.10 

239 

-0.017 

X 

4 

e 

0.97 

2.44 

0.11 

157 

0.0 

X 

3 

f 

0.99 

2.50 

0.11 

141 

0.021 

X 

3 

g 

1.02 

2.57 

0.11 

114 

0.047 

X 

3 

h 

1.06 

2.66 

0.11 

103 

0.084 

X 

4.5 

i 

1.10 

2.77 

0.11 

144 

0.129 

-0.0125 

7 

a 

0.76 

1.98 

0.18 

59 

-0.129 

-0.390 

7 

b 

0.83 

2.14 

0.18 

64 

-0.066 

X 

X 

6 

c 

0.89 

2.30 

0.19 

86 

0.0 

5.5 

d 

0.95 

2.46 

0.19 

99 

0.074 

X 

4 

e 

1.00 

2.59 

0.20 

59 

0.135 

X 

6 

f 

1.04 

2.68 

0.21 

64 

0.183 

-0.390 

8 


Table 2 - Control Derivatives and Transfer Function Zeros for Meeker-Hall Group AB. 

(Data from Meeker and Hall 1967) 
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BB-l 

BB-2 

BB-3 

(Odr (rad/sec) 

2.49 

2.49 

2.60 

u 

0.10 

0.10 

0.08 

l<t>/Pldr 

5.22 

4.82 

4.90 

Z«t>/p) dr (deg) 

42.0 

46.0 

47.0 

^roll (sec) 

0.4 

0.37 

0.4 

^spiral (sec) 

997.0 

987 

166 


Table 3 - Modal Data for Meeker-Hall Group BB. 
(Data from Meeker and Hall 1967) 


Sub- 

Group 


BB-1 


BB-2 


BB-3 


Config 

Cty/COdr 

coq) 

(rad/sec) 

Co 

L'5 

(deg/sec 2 )/in 

N’t/L's 

N 'p 

(sec* ) 

Pilot 

Rating 

a 

0.48 



115 

—iIiL’Hi 

0.0473 

7 

b 

0.73 

1.81 

0.06 

160 

-0.054 

si 

4 

c 

0.73 

1.81 

0.06 

153 

-0.054 

si 

6.5 

d 

0.76 

1.89 

0.06 

190 

-0.047 

si 

7 

e 

0.85 

2.11 

0.07 

71 

-0.027 

si 

4 

f 

0.96 

2.38 

0.07 

116 

0.0 

si 

5 

g 

0.96 

2.38 

0.07 

152 

0.0 

1 

6 

h 

1.14 

2.83 

0.07 

72 

0.052 

0.0473 

8 

a 

0.66 

1.65 

0.08 

88 

-0.066 

0.01477 

6.5 

b 

0.81 

2.02 

0.09 

125 

-0.034 

l 

i 

3 

c 

0.94 

2.34 

0.09 

190 

0.0 

2 

d 

1.06 

2.65 

0.1 

95 

0.037 

si 

4 

e 

1.18 

2.93 

0.1 

76 

0.074 

0.01477 

7 

a 

0.23 

0.59 

0.15 

91 

-0.111 

-0.0989 

7.5 

b 

0.66 

1.72 

0.14 

185 

-0.047 

si 

4 

c 

0.85 

2.20 

0.16 

207 

0.0 

si 

3 

d 

0.97 

2.52 

0.18 

160 

0.037 

si 

2 

e 

1.08 

2.80 

0.19 

95 

0.074 

si 

4 

f 

1.19 

3.09 

0.21 

116 

0.116 

si 

3.5 

g 

1.35 

3.50 

0.23 

54 

0.183 

-0.0989 

5.5 


Table 4 - Control Derivatives and Transfer Function Zeros for Meeker-Hall Group BB. 

(Data from Meeker and Hall 1967) 
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CB-1 

CB-2 

CB-3 

Ci>dr (rad/sec) 

2.47 

2.48 

2.44 

Cdr 

0.10 

0.10 

0.08 

WPldr 

12.0 

14.0 

14.7 

Z(4>/p) dr (deg) 

41.0 

39.0 

47.0 

Troll (sec) 

0.4 

0.5 

0.4 

Tspiral (sec) 

-250 

111 

-104 


Table 5 - Modal Data for Meeker-Hall Group CB. 
(Data from Meeker and Hall 1967) 


Sub- 

Group 


CB-1 


CB-2 


CB-3 


Config. 

Oty/COdr 

(rad/sec) 

** 

L'8 

(deg/sec 2 )/in 

N'tfL's 

N \ 
(sec ) 

Pilot 

Rating 

a 

~ 

-1.54* 

1.49* 

70 

-0.066 

0.0351 

9 

b 

0.42 

1.03 

0.04 

153 

-0.034 

4 

4 

7 

c 

0.87 

2.15 

0.03 

110 

0.0 

6.5 

d 

1.01 

2.49 

0.03 

151 

0.015 

4 

7 

e 

1.18 

2.91 

0.03 

112 

0.037 

0.0351 

9 

a 

— 

-1.68* 

1.80* 

156 

-0.066 

0.0180 

10 

b 

0.28 

0.70 

0.12 

107 

-0.034 

4 

4 

4.5 

c 

0.28 

0.70 

0.12 

115 

-0.034 

7 

d 

0.28 

0.70 

0.12 

115 

-0.034 

4 

5 

e 

0.83 

2.06 

0.11 

103 

0.0 

4 

5.5 

f 

1.16 

2.88 

0.14 

110 

0.037 

0.0180 

8 

a 

— 

-0.92* 

0.96* 

146 

-0.034 

-0.0108 

7.5 

b 

0.72 

1.75 

0.14 

150 

0.0 

4 

6 

c 

1.11 

2.71 

0.19 

110 

0.037 

4 

5.5 

d 

1.29 

3.14 

0.22 

107 

0.059 

4 

8 

e 

1.29 

3.14 

0.22 

146 

0.059 

4 

8 

f 

1.49 

3.62 

0.25 

64 

0.087 

-0.0108 

8 


* denotes real zeros 


Table 6 - Control Derivatives and Transfer Function Zeros for Meeker-Hall Group CB. 

(Data from Meeker and Hall 1967) 
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AB-1 

AB-2 

AB-3 

g/V 0 (sec-1) 

0.0547 

0.0547 

0.0547 

L) 3 (sec-2) 

-11.20 

-13.34 

-10.58 

L' p (sec-1) 

-2.66 

-2.49 

-2.20 

L' r (sec-1) 

0.901 

0.730 

1.195 

N'p (sec-2) 

5.935 

6.012 

5.289 

N' r (sec-1) 

-0.2240 

-0.3510 

-0.7061 

Yp’ (sec-1) 

-.1594 

-0.1616 

-0.1559 

Y p~ l 

Yn + a-o 

-1.014 

-1.014 

-1.016 

0.0207 

0.0238 

0.0251 

Yr- 1 

-0.9980 

-0.9966 

-0.9928 

J 

| BB-1 

BB-2 

BB-3 


g/V 0 (sec-1) 

0.0547 

0.0547 

0.0547 

L'p (sec-2) 

-45.12 

-42.14 

-40.62 

L'p (sec-1) 

-2.68 

- 2.74 

-2.21 

L’ r (sec-1) 

1.281 

2.058 

4.422 

N'p (sec-2) 

5.738 

5.544 

4.857 

N' r (sec-1) 

-0.169 

-0.278 

-0.575 

Yp’ (sec-1) 

-0.156 

-0.159 

-0.147 

Y r l 

-1.015 

-1.016 

-1.014 

Y p + a 0 

0.0219 

0.0207 

0.0243 

Yr- 1 

-0.998 

-0.997 

-0.994 



CB-1 

CB-2 

CB-3 

g/V 0 (sec-1) 

0.0547 

0.0547 

6.0547 

L'p (sec-2) 

-106.4 

-110.9 

-116.2 

L'p (sec-1) 

-2.864 

-2.237 

-2.414 

L' r (sec-1) 

-0.0626 

7.768 

13.96 

N'p (sec-2) 

4.689 

4.265 

3.062 

N' r (sec-1) 

0.0134 

-0.3190 

-0.3451 

Yp’ (sec-1) 

-0.1420 

-0.1496 

-0.1283 

Y p ~ 1 

-1.015 

-1.016 

-1.015 

Yp + a o 

0.0225 

0.0198 

0.0225 

Y r - 1 

-0.998 

-0.997 

-0.997 


Table 7 - Stability Derivatives Meeker-Hall Groups AB, BB, and CB. 
(Data from Meeker and Hall 1967) 
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l<t>/Pldr 

0 

^(<KP)dr (degrees) 
60 

120 

0 

✓ 

X 

X 

1.5 

✓ 

✓ 

✓ 

3.0 

✓ 

✓ 

✓ 


V denotes tested configuration 
X denotes not tested 


Table 8 - Costigan-Calico Test Matrix. 
(Data from Costigan and Calico 1989) 


State Vector 


Roll Mode 


Spiral Mode 


Dutch Roll 
Mode 


Eigenvalue 


-4.0 


T 

P 

r 


-0.025 


( 3.0 , 0.4 )* 


Eigenvector 


0 

1 

0 


* (•,•) denotes (frequency (rad/sec), damping) for eigenvalues 
x denotes elements not weighted in the cost function 


0 

x 

X 

1 


Table 9 - Costigan-Calico Design Parameters. 
(Data from Costigan and Calico 1989) 


TGT 

X 

X 

^dr 


Kp/Pldr 

Yaw Pointing Task 
^(<t>/P)dr (degrees) 

0 120 

Bank Angle Tracking Task 
^(<KP)dr (degrees) 

0 120 

0 

6.2 

- 

4.2 

.. 

1.5 

5.2 

5.7 

4.0 

4.5 

3.0 

4.2 

6.0 

3.8 

4.7 


Table 10 - Costigan-Calico Tracking Task Average Cooper-Harper Ratings. 
(Data from Costigan and Calico 1989) 
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